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Transverse quasilinear relaxation of the cyclotron-Cherenkov instability 
in the inhomogeneous magnetic field of pulsar magnetospheres is considered. 
We find quasilinear states in which the kinetic cyclotron-Cherenkov instability 
of a beam propagating through strongly magnetized pair plasma is saturated 
by the force arising in the inhomogeneous field due to the conservation of the 
adiabatic invariant. The resulting wave intensities generally have nonpower 
law frequency dependence, but in a broad frequency range can be well approx- 
imated by the power law with the spectral index —2. The emergent spectra 
and fluxes are consistent with the one observed from pulsars. 



I. INTRODUCTION 

In this paper we consider quasilinear relaxation in inhomogeneous magnetic field of a 
highly relativistic beam propagating along the strong magnetic field through a pair plasma. 
This describes the physical conditions on the open field lines of pulsar magnetospheres 
(e.g. [IH). The possibility of the cyclotron-Cherenkov instability of the beam in the pulsar 
magnetosphere has been suggested by [@] and and developed later by ||5|, 0. The 
cyclotron-Cherenkov instability develops at the anomalous Doppler resonance 

uj{\l) - hm\ - s— = for s < (1) 

7 

where co'(k) is the frequency of the normal mode, k is a wave vector, v is the velocity of 
the resonant particle, ujb = \q\B/mc is the nonrelativistic gyrofrequency, 7 is the Lorentz 
factor in the pulsar frame, q is the charge of the resonant particle, m is its mass and c is the 
speed of light. It has been shown (e.g. [^), that cyclotron-Cherenkov instability can explain 
a broad variety of the observed pulsar phenomena. 
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Close to the stellar surface, where the initial beam is produced and accelerated, the 
particles quickly reach their ground gyrational state due to the sychrotron emission in a 
superstrong magnetic field, so that their distribution becomes virtually one dimensional [|I| . 
In the outer parts of magnetosphere it becomes possible to satisfy the anomalous Doppler 
resonance - the cyclotron- Cherenkov instability develops bringing about the diffusion of 
particles in transverse moments. The relevant saturation mechanism then determines the 
final spectrum (which can be later modified be the absorption processes). 

The nonlinear saturation of the cyclotron-Cherenkov instability due to the diffusion of 
the resonant particles has been previously considered by several authors. Kawamura & 
Suzuki neglected the possible stabilizing effects of the radiation reaction force due to the 
cyclotron emission at the normal Doppler resonance and the force arizing in the inhomoge- 
neous magnetic field due to the conservation of the adiabatic invariant. These forces result 
in a saturation of the quasilinear diffusion. 

Lominadze et al. were the first to notice the importance of the radiation reaction force 
due to the emission at the normal Doppler resonance on the saturation of the quasilinear 
diffusion. Unfortunately, Lominadze et al. used an expression for the cyclotron damping 
rates which is applicable only for the nonrelativistic transverse motions, when 7?/' {psi is 
the pitch angle) is much less than unity. In the pulsar magnetospheres the development of 
the cyclotron-Cherenkov instability results in a diffusion of particles in transverse moments, 
quickly increasing there transverse energy to relativistic values. 

In a review paper Lominadze et al. took a correct account of the radiation reaction 
force due to the emission at the normal Doppler resonance and pointed out the importance 
of the the force arizing in the inhomogeneous magnetic field due to the conservation of 
the adiabatic invariant ( G force Eq. (p!OD). When considering the deceleration of the beam 
Lominadze et al. [0 has incorrectly neglected the radiation reaction force due to the emission 
at the anomalous Doppler resonance in comparison with the radiation reaction force due to 
the emission at the normal Doppler resonance. 

In this work we reconsider and extend the treatment of the quasilinear stage of the 

2 



cyclotron-Cherenkov instability. We have found a state, in which the particles are constantly 
slowing down their parallel motion, mainly due to the component along magnetic field of 
the radiation reaction force of emission at the anomalous Doppler resonance. At the same 
time they keep the pitch angle almost constant due to the balance of the force G and the 
component perpendicular to the magnetic field of the radiation reaction force of emission 
at the anomalous Doppler resonance. We calculate the distribution function and the wave 
intensities for such quasilinear state. 

In the process of the quasilinear diffusion the initial beam looses a large fraction of its 
initial energy ^ 10%, which is enough to explain the typical luminosities of pulsars. Though 
the quasisteady wave intensities are not strictly power laws (see Eq. (|38|)), they can well 
approximated by a power law with a spectral index -F(i^) oc —2 (-F(z^) is the spectral flux 
density in Janskys) which is very close to the observed mean spectral index of —1.6 [§. The 
predicted spectra also show a turn off at the low frequencies v < 300MHz and a flattering 
of spectrum at large frequencies u > IGHz which may be related to the possible turn-up in 
the flux densities at mm- wavelengths 

II. QUASILINEAR DIFFUSION 

A particle moving in a dielectric medium in magnetic fleld with the velocity larger than 
the velocity of light in a medium is emitting electromagnetic waves at the anomalous Doppler 
resonance ( s < in Eq. (|l|)) and at the normal Doppler resonance ( s > in Eq. (^). The 
radiation reaction due to the emission at the normal Doppler resonance slows the particle's 
motion along magnetic fleld and decreases its transverse momentum. The radiation reaction 
due to the emission at the anomalous Doppler resonance increases its transverse momentum 
and also slows the particle's motion along magnetic fleld |T^. As the particle propagates 
into the region of lower magnetic fleld, the force G decreases the its transverse momentum 
and increases the parallel momentum. The stationary state in transverse moments may be 
reach when the actions of the G force and radiation reaction due to the emission at the 
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normal Doppler resonance is balanced by the radiation reaction due to the emission at the 
anomalous Doppler resonance. 

The quasisteady stage may also be considered in terms of a detailed balance of for the 
particle transitions between the Landau levels. The quasisteady stage is reached when the 
number of induced transitions up in Landau levels due to the emission at the anomalous 
Doppler resonance is balanced by the number of the spontaneous transitions down in Landau 
levels due to the emission at the anomalous Doppler resonance. 

The equations describing the quasilinear diffusion in the magnetic field are 
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where E'^(k)d\<./ (27r)^ is the energy density of the waves in the unit element range of k-space. 
In the Table | we give the dimensions of the main used quantities. 

In Eq. (|^) we neglected the spontaneous emission processes at the anomalous Doppler 
resonance and the induced emission processes at the normal Doppler resonance. The net 
effect of the spontaneous emission at the normal Doppler resonance is treated as a damping 



force acting on each particle in the Boltzman-type left hand side of equation (^. To be 
exact, we should have treated the effects of spontaneous emission at the normal Doppler 
resonance as stochastic terms in the Fokker-Plank-type terms on the right hand side of 
equation (|^). But the fact that the emission at the normal Doppler resonance occurs on 
very high frequencies at which the presence of a medium is unimportant in the dispersion 
relation of the waves and can be neglected allows one to integrate the corresponding terms 
over angles and sum over harmonics to obtain a classical synchrotron radiation damping 
force, that can be treated using the Boltzman approach. Thus, the total time derivative of 
the distribution function is 



dfip) dfip) , ^9/(p) ^ d_ 



G + F + ^(vxBo))/(p) 



(9) 



dt dt dr dp 

where G is the force due to the conservation of the adiabatic invariant 
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G|| = -/37^^ Gx = _/37^, /3 = ^ (10) 

Here Rb ^ 10^ cm is the radius of curvature and F is the radiation damping force due to 
the spontaneous synchrotron emission at the normal Doppler resonance: 

F||=-«7V, F^ = -a^(l + 7V), « = ^ (11) 
From (|10]) and (0) we find that 

^ = ^7, ^ = ;^^^' for7V^>l (12) 

where ri = c/uob is a Larmor radius and re = /{mc^) = 2.8 x 10~^^cm is a classical radius 
of an electron. 

The dimensionless ration in (|T2|) is 



5 X IQ-^ RB,9Rg^ (13) 



' L 

Rb,9 = -Rfi/lO^cm is the radius of curvature in units of 10^ cm, Rg = /2/lO^cm is the distance 
from the neutron star surface in units of 10^ cm. 



Using ([I3D we find that for tlie primary particles witli 7 ^ 10''' 
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Tlien the total derivative (|^) may be written as 

df{p) <9/(p) (9/(p) 1 ^ / • ,^ 1 / 2p N 



(15) 



dt dt dr p sin -0 dijj p"^ dp 

We are interested in the quasilinear diffusion of the particles due to the resonant interac- 
tion with the waves at the anomalous Doppler effect. We expand the transition currents (H) 
in small v± keeping only s = — 1 terms: V(— l,p, k) = v±/2 {l,ia,0). Then for the waves 
propagating along magnetic field e(k) = (1, 0, 0) we find 



w{±l, p, k) = -^[J^^i^O^) ~ s^b/i - hn) 



(16) 



where we took into account that -R£;(k) 1/2. 

We now can find the diffusion coefficients in the approximation of a one dimensional 
spectrum of the waves. 
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We first note that we can simplify the change in the pitch angle (|^) in the limit ip"^ 5 



and 1/7^ < 5 
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We then find 
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where 



El = hu{k)n{k) 
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is energy density per unit of one-dimensional wave vector and we assumed that ci;(k) is an 
isotropic function of k. 

We next solve the partial differential equation describing the evolution of the distribution 
function by sucsessive approximations. We first expand equation (|^) in small ip assuming 
that 8/ dip ^1/?/'. We also neglect the convection term assuming that the characteristic 
time for the development of the quasilinear diffusion is much smaller that the dynamic time 
of the plasma flow. Then we assume that it is possible to separate the distribution function 
into the parts depending on the ijj and p: 
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with 
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In the lowest order in we obtain an equation: 
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The next order in ip gives 
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By integrating (^) over ip with a weight ip we find the equation for the parallel distri 
bution function: 



(26) 



where 

The term containing A describes the slowing of the particles due to the radiation reaction 
force and the term containing D describes the slowing of the particles due to the quasilinear 
diffusion, or, equivalently, due to the radiation reaction force of the anomalous Doppler 
resonance. To estimate the relative importance of these terms we consider a ratio 

D^c-^- ^ ^ ^^^^ 

Neglecting the second term on the left hand side of (^) we find 

^ - 11 {pD'"^''='Einp)) == (29) 

If the cyclotron quasilinear diffusion has time to fully develop and reach a steady state, 
then the distribution function of the resonant particles is 



fip) « (30) 



Next we turn to the equation describing the temporal evolution of the wave intensity 
. Neglecting the spontaneous emission term and the wave convection we find 



dt 

where 
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and we introduced 

filh'd^ = fip)p'dp (33) 

We will estimate this growth rate for the emission along the external magnetic field for 
distribution (P^D, (^). Neglecting d/dp and assuming that ip^ <^ 26 (so that most of the 



particles are moving with the superluminal velocity) we find for s = — 1 
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(It is important to note that in the hmit ip"^ -C 25 the growth rate does not depend on the 
scatter in pitch angles). 

Equations (|29|) and (^ID may be combined to a quasihnear expression 

Which after integration gives 

Neglecting the initial density of particles in the region of quasilinear relaxation and using 
Eqs. ( ^0|) and (|36|) we can find a distribution function and the asymptotic spectral shape: 

1/2 

Uog(7maa;/7) log(7max/7mm) 



/(7) = ,..L.. ^ 1 (37) 



^ 27Crerl \log{'J^ax/'lmin) J 2'KUj'^r \ logi^f max hmin) ) 

It is noteworthy, that a simple power law distribution for the spectral intensity and 
distribution function cannot satisfy both Eqs. (|30|) and (|36D . The particle distributions 
function and the energy spectrum of the waves are displayed in Figs. |I| and 0. 
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FIG. 1. Asymptotic distribution functions in j — p± and 'j — ip spaces in arbitrary units 
for jmax = 10''. The spike at the 7 = jmax is an artifact of the initial distribution function 
fi'y)^ — '^(7 ~ Iraax)- The divergence at 7 = ^max is weak (logarithmic) and would be removed if 
the more realistic nitial distribution function was used. 
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FIG. 2. Asymptotic one dimensional energy density in the waves in the 7-space (arbitrary 

units), and the predicted observed flux in Janskys. 



We can now estimate the flux per unit frequency: 



log (7^ 



(39) 
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characteristic pitch angle 

^o = s( ^^^] { ^"g^^— \ ^ lo-^ (40) 

\ rl ) \\oginmaxhmin) j 

(which remarkably stays almost constant for a broad range of particles' energies and also 
for different values of 7mm), and the total energy density in the waves 

Etot = / F{u)du ^ -, (41) 

•^''min 4Vvrrer|.log / (7max/7mm) 

This total energy can be compared with the kinetic energy density of the beam: 

^ / (42) 

It means that some considerable fraction of the beam energy can be transformed into waves. 

We can also estimate the energy flux (|39|) at the Earth. Assuming that distance to the 
pulsar is 1 kpc, we find 

F°^'iv) ^ 300 Jy ( ^-—] ^ (43) 

^ ^ ^ V400MHzy ^ ^ 

With time, the value of 7mm decreases as the particles are slowed down by the radiation 
reactions force. Since at the given radius, the particles with lower energies resonate with 
waves having larger frequencies, more energy will be transported to higher frequencies hard- 
ening the spectrum. The lower frequency cutoff is determined by the initial energy of the 
beam. No energy is transported to frequencies lower than 

UJmin = —7 (44) 

This simple picture, of course, will be modified due to the propagation of the flow in the 
inhomogeneous magnetic field of pulsar magnetosphere. 



III. CONCLUSION 



In this work we investigated the new saturation mechanism for the cyclotron-Cherenkov 
instability of a beam in a inhomogeneous magnetic field. We showed that for the typical 
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parameters of the pulsar magnetosphere it is possible to reach quasisteady state, in which 
the transverse motion of particles is determined by the balance of a radiation reaction force 
due to the emission at anomalous Doppler effect and the force arising in the inhomogeneous 
magnetic field due to the conservation of adiabatic invariant. The resulting wave intensities 
are sufficient to explain the observed fluxes from radio pulsars. 
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